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We compare ferromagnetic resonance measurements of Permalloy Ni8oFe20 (Py) films sputtered 
onto Cu(lll) films with and without a graphene (Gr) interlayer grown by chemical vapor deposition 
before Py deposition. A two-angle sputtering method ensured that neither Gr nor Py was degraded 
by the sample preparation process. We find the expected damping enhancement from spin pumping 
for the Py/Cu case and no detectable enhancement for the Py/Gr/Cu case. Since damping is 
sensitive to effects other than spin pumping, we used magnetometry to verify that differences in 
Py magnetostatic properties are not responsible for the difference in damping. We attribute the 
suppression of spin pumping in Py/Gr/Cu to the large contact resistance of the Gr/Cu interface. 


I. INTRODUCTION 

The complex interactions between magnetization dy¬ 
namics in a ferromagnet and the flow of charge and spin 
currents is a topic that combines rich physics and rel¬ 
evance to technological applications. These interactions 
are often divided into two types that are reciprocal mani¬ 
festations of the same microscopic process [1]. Spin trans¬ 
fer refers to the torque exerted on the magnetization by 
a flow of spins (with or without an associated flow of 
charge). Spin pumping refers to the flow of spins gener¬ 
ated by a precessing magnetization. The first effect en¬ 
ables control of magnetization through an applied current 

[ 2 ] , rather than an applied magnetic field, and is crucial 
to realizing a scalable magnetic random access memory 

[3] . The second effect is most often encountered as an en¬ 
hanced damping of magnetization dynamics when a thin 
ferromagnetic (FM) film is in contact with a nonmagnetic 
(NM) material that absorbs the spin current [4, 5]. This 
loss of angular momentum adds to the intrinsic dissipa¬ 
tion of the FM material, causing an enhancement of the 
Gilbert damping parameter a. The effect on damping 
can be quite large 2x) for ultrathin (< 10nm) FM 
films and has been observed for many FM/NM combina¬ 
tions [5-11]. 

Brataas et al. [1] have reviewed the phenomenology 
and the quantitative theory of spin transfer and spin 
pumping. When considering the effect of spin pumping 
on FM damping, it is useful to identify two limiting cases. 
If the FM is in contact with a material that strongly scat¬ 
ters or absorbs the pumped spin current, a is increased 
substantially. In contrast, if the NM does not act as a 
spin sink, then pumped spins will accumulate in the NM 
and drive a diffusive spin current back toward the FM 
that cancels the spin pumping current in steady state. 
In this “spin battery” limit, there is no effect on a. For 
a general case, the change in damping for a FM film of 
thickness dpM is determined by an effective spin conduc¬ 
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tance per unit area, Geft (units of D ^ • m ^), according 
to 


Ac^sp 


l7l^^ 

2e^MsdFM 


( 1 ) 


where 7 is the gyromagnetic ratio, e is the electron 
charge, H is the reduced Planck constant, and Mg is the 
saturation magnetization of the FM. The conductance 
GeS includes contributions from the FM/NM interface 
and from spin diffusion in the NM material (see Ap¬ 
pendix A). 

Graphene (Gr) is expected to be a useful material for 
spintronics because its high carrier mobility and low spin 
scattering should allow propagation of spin currents over 
long distances (^ 1 pm) at room temperature. Spin co¬ 
herence lengths of about 1 pm have been demonstrated 
by injecting spin-polarized charge current from FM con¬ 
tacts into a Gr channel [12]. The main challenge facing 
these devices is the limited efficiency for injecting spin- 
polarized carriers from the metallic FM into a channel 
having much lower electrical conductivity. Even with 
tunnel barriers between the FM and the Gr designed to 
minimize this effect, the efficiency is only about 30% ]12]. 
Since spin pumping involves only spins and not charges 
flowing between the FM and NM materials, it provides a 
way to inject spin current from FM metals into Gr that 
is not limited by the mismatch in electrical conductivity. 

This paper compares ferromagnetic resonance (FMR) 
measurements of Permalloy Ni 8 oFe 2 o (Py) films with 
nominal thicknesses between 8 nm and 45 nm sputtered 
onto Gu(lll) films with and without a Gr interlayer 
grown by chemical vapor deposition (GVD) before Py de¬ 
position. The Py films were deposited using a two-step 
process (see section III A) that avoids damage to the Gr 
and yields high quality Py. The FMR results show the ex¬ 
pected damping enhancement for the Py/Gu case but no 
enhancement for the Py/Gr/Cu case. This unexpectedly 
strong effect on spin tranport by a monolayer of material 
with weak intrinsic spin scattering requires careful val¬ 
idation. Since FM damping is sensitive to effects other 
than spin pumping, much of the paper is concerned with 
ensuring that neither Gr nor Py is significantly degraded 
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by the sample preparation process. In particular, we use 
a novel two-angle sputtering geometry, we verify the Gr 
quality using Raman spectroscopy, we use magnetometry 
to measure the moment and coercivity of each sample, 
and we examine the Py morphology using transmission 
electron microscopy. 

Previous studies have shown evidence for spin pumping 
from Py into Gr [13, 14], but the use of only a single 
thickness of Py and the lack of magnetometry or other 
characterization of the Py itself limits the conclusions 
that can be drawn from the observed effects. We discuss 
these limitations in detail, drawing on our results and on 
other recently published work [15]. Finally, we attribute 
the suppression of spin pumping in Py/Gr/Gu to the 
large contact resistance of the Gr/Gu interface [16[. 


0.014 

0.012 

D) 

g 

CL 

TO 0.010 

Q 

a 

0.008 


0.006 

0 10 20 30 40 50 

4 m (nm) 


i 


■ Py/Cu 

■ Py/Gr/Cu 

a 


[ ^ 


1 - 

H 1 



II. FMR MEASUREMENT OF SPIN PUMPING 

For FMR measurements, samples were placed on a 
coplanar waveguide with a center conductor width of 
100 /im, and the scattering parameter S 21 was measured 
using a vector network analyzer. The static magnetic 
field H was applied perpendicular to the sample plane 
and was large enough to saturate the magnetization, 
suppressing the 2 -magnon scattering contribution to the 
FMR linewidth [17[. For each excitation frequency /, 
fitting the real and imaginary parts of S 2 i{H) yields 
the resonance field Fires (/) and the field-swept linewidth 
AH{f). Details of the measurement and data analysis 
techniques are described in [18]. 

The FMR resonance field and linewidth as a function 
of excitation frequency are given by 

Plres(/) = ^+Meff (2) 

AH{f) = ^ + AHo, (3) 

ItImo 

where po is the magnetic constant, MgS = Mg — 
is the effective magnetization including a perpendicular 
anisotropy field and AHq is the broadening due to 
inhomogeneity in the local resonance field across the sam¬ 
ple [18]. A linear fit to Eq. 2 gives values for 7 and Mgg, 
and a linear fit to Eq. 3 gives values for a and AHq 
[18]. Results for all four fitted parameters are given in 
Appendix B. 

Figure 1 shows a as a function of dpM for two sets of 
Py samples. The Py/Gu samples clearly show an increase 
in a as dpM decreases, as seen previously [ 8 ] and as ex¬ 
pected from Eq. 1, while the Py/Gr/Cu samples show no 
significant change in a. As shown in the next section, the 
Py in the two sets of samples has very similar magnetic 
and morphological properties, so this result is strong ev¬ 
idence for suppression of spin pumping between Py and 
Gu due to the Gr interlayer. 


FIG. 1. Damping vs. thickness for Py deposited on bare Cu 
(Py/Cu) and on Cu covered by CVD Gr (Py/Gr/Cu). Ver¬ 
tical error bars show the uncertainty from the linear fits to 
Eq. 3. The larger uncertainty in a for the thicker Gr sam¬ 
ples is due to a larger scatter in AH{f) for these samples. 
Horizontal error bars reflect an estimated 10% uncertainty in 
deposition rate. 


III. SAMPLE PREPARATION AND 
CHARACTERIZATION 

A. Motivation for Two-angle Sputtering 

Damage to Gr during magnetron sputtering of Ti and 
A1 has been studied using Gr flakes exfoliated onto Si sub¬ 
strates [19]. This work considered several types of bom¬ 
bardment present in sputtering (electrons, photons, de¬ 
positing atoms, and inert gas atoms) and concluded that 
damage is primarily due to inert gas atoms, which can be 
reflected from the target as neutral atoms with sufficient 
energy to displace C atoms from Gr. Detailed simulations 
[20] of sputtering a Cu target in Ar gas show, for a typical 
discharge voltage of 400 V, that Cu atoms leave the tar¬ 
get with an average energy of about 10 eV, while Ar neu¬ 
trals are reflected with an average energy of about 45 eV. 
For comparison, the displacement threshold for Ar strik¬ 
ing Gr at normal incidence is 33 eV [21, 22[. Increasing 
the gas pressure or target-to-substrate distance reduces 
the flux of both Ar neutrals and depositing atoms, but 
since the more energetic Ar atoms have a longer mean 
free path than the depositing atoms [ 20 [, the deposition 
rate is strongly suppressed before the Gr is protected 
from Ar impacts. Since the displacement threshold in¬ 
creases away from normal incidence [ 21 , 22 ], damage can 
be avoided by orienting the substrate perpendicular to 
the target, as seen in Fig. 2 with 9 = 90°. The tradeoff 
between deposition rate and Gr damage for orientations 
near perpendicular was quantified for Ti and A1 in [19]. 

For many magnetic materials, including Py, bombard¬ 
ment of the substrate is important in promoting growth 
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of thin films that are homogeneous and have properties 
such as Mg and a that are close to bulk values. Thus a 
substrate orientation perpendicular to the target is unde¬ 
sirable in terms of film quality as well as deposition rate. 
For this reason, we chose to deposit Py at two separate 
orientations, as shown in Fig. 2. The first deposition, 
5 nm at 0 = 90°, covers the Gr with a Py film of low 
quality but does not damage the Gr, and this Py layer 
protects the Gr from damage during the second deposi¬ 
tion, at 0 = 30°, that yields high quality Py. The two 
depositions were done without breaking vacuum by use 
of an angled sample holder (see Methods section). Im¬ 
portantly, transmission electron microscopy (see Section 
IIIG and Fig. 5) shows the final Py film to be a sin¬ 
gle, homogeneous layer rather than two distinct layers. 
We attribute the lack of a morphological boundary be¬ 
tween the two layers to bombardment during the second 
deposition. The sum of the first and second deposition 
thicknesses is the nominal Py thickness, do- 



FIG. 2. Schematic of deposition at two angles 0 relative to 
substrate normal direction (not to scale). The Py target is 
shown in different positions for clarity, but in practice 0 was 
changed by tilting the substrate. 


B. Gr Properties 

We used Raman spectroscopy to measure Gr damage 
as follows. The Raman spectrum of pristine Gr has two 
characteristic peaks, named G and 2D, and defects allow 
an additional scattering mechanism that adds a feature 
named the D peak [23]. The ratio of the peak intensities 
for the D and G features. Id/Iqj can be used to estimate 
the number of point defects per unit area [24]. 

For our as-grown Gr on Cu, we found a baseline value 
of /d//g = 0.11 ± 0.03, which implies one defect ev¬ 
ery (1300 ± 350) nm^. We deposited 5 nm to 35 nm of 


Py at 9 = 90° on several samples and found no change 
in /d//g) demonstrating that sputtering in this geom¬ 
etry does not damage Gr. The fact that Gr Raman 
peaks could be measured through these films indicates 
the Py deposited at 6 = 90° was highly transparent, con¬ 
sistent with the fact that these samples still showed some 
Gu color after deposition. In contrast, Py deposited at 
6 = 30° did not show any Gr Raman peaks for thickness 
> 2 nm and had a color similar to that of high quality 
Ni films. A short deposition at 6 = 30°, lasting w 2s 
and giving Ri 0.2 nm, allowed us to measure Gr Raman 
peaks with Id/Ig = 1-7, corresponding to one defect ev¬ 
ery 100 nm^. This confirms that Gr is readily damaged 
by sputtering in this geometry, with a damage rate of 
roughly 0.005nm“^ • s“^. In fact, this rate is close to the 
flux of Ar neutrals expected from the simulations in [20] 
(see Appendix C). 

We estimated the potential damage to Gr covered by 
Py using the Stopping and Range of Ions in Matter 
(SRIM) Monte Carlo simulation [25]. Since the substrate 
in our case is located 23 cm from the target, only Ar 
neutrals reflected nearly perpendicular to the target will 
reach the substrate, with a maximum energy of « 70 eV 
for a 430 V discharge [20]. For 10^ Ar ions at 70 eV 
and normal incidence, SRIM predicts no transmission for 
5 nm of Py with its nominal density of 8.7g/cm^. Ad¬ 
justing parameters to reflect the fact that Py deposited at 
9 = 90° is different from nominal, e.g., using a thickness 
of 2 nm and nominal density, or a thickness of 5 nm and 
a density of 3g/cm^, yielded a transmission probability 
of less than 1 in 10^. Combining this with the Ar neutral 
flux estimated in Appendix C shows that damage due to 
Ar penetration through the 9 = 90° Py film is expected 
to be negligible. 


C. Py properties 

In addition to FMR measurements, we used vibrating 
sample magnetometry (VSM) to measure Py magneto¬ 
static properties of the same samples presented in Fig. 1. 
We also used transmission electron microscopy (TEM) to 
characterize Py morphology for separate, similarly pre¬ 
pared samples. 

For a homogeneous FM film with magnetization Mg 
and thickness dg, the moment per unit area is /ioMsdg. 
The total magnetic moment measured by VSM, divided 
by the area of each sample, is plotted vs. nominal Py 
thickness in Fig. 3. From the slopes of the linear fits 
in Fig 3, we find /rgMs = (1.1 ± 0.1) T for both Py/Cu 
and Py/Gr/Cu samples. The nonzero intercepts of the 
fits imply the actual thickness of the films is « 3nm 
less than the nominal thickness. This is not surprising, 
since the 5 nm of Py deposited at 9 = 90° is porous (see 
Section IIIB) and some of the Py deposited at 9 = 30° 
will fill the pores instead of adding to the overall film 
thickness. For the rest of this paper, we use the actual 
FM thickness dpM obtained by subtracting the intercept 
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values from do. 

Given the actual FM thickness values, we can use 
Meff measured by FMR to determine the perpendicu¬ 
lar anisotropy of the Py as follows. Figure 4 shows 
plots of Meff vs. f/dpM for the two-angle Py/Cu and 
Py/Gr/Cu samples, and for several samples deposited 
only at 6* = 30°. The linear dependence for all samples is 
consistent with a constant interfacial anisotropy energy 
for which Mgff is given by 


Meff = Ms — 


2Kint 

/iO-^sdpM 


Thus the intercept of each linear ht in Fig. 4 gives Ms 
and the slope gives Kint- We find fXoMg = (1.01 ± 0.02) T 
for both Py/Gu and Py/Gr/Gu samples, consistent with 
the value obtained from VSM data. The values of Ki^t, 
given in Table 1, are quite similar for a given deposi¬ 
tion method, but are noticeably larger for two-angle films 
than for single-angle films. Furthermore, our values are 
larger than the value of Ki^t « 1 x 10“‘*J/m^ reported 
for Py/Gu [26] and Py/air [27] interfaces. These differ¬ 
ences are likely due to variations is how Py nucleates and 
grows on different substrates, for example polycrystalline 
Cu compared to our Cu(lll), in addition to the different 
deposition geometries. 

The coercive field measured by VSM is 0.5 mT for 
Py/Cu and 0.25 mT for Py/Gr/Cu. Although both val¬ 
ues are somewhat larger than the w 0.1 mT typically 
found for high-quality Py, similar values have been re¬ 
ported for Py deposited on Cu(lll) single crystals [28]. 
Finally, a 40 nm Py/Cu film deposited at a single angle 
oi 9 = 30° showed the same value of a (w 0.008) as the 
two-angle Py/Cu film with 5 nm at 0 = 90° and 40 nm 
at 0 = 30°, indicating that the intrinsic damping of the 
Permalloy is not affected by two-angle deposition. 

Overall, these results show that whatever effects Gr 
may have on the film growth process, for our two-angle 
method it has little or no effect on the Py magnetic prop¬ 
erties. 

We prepared samples for cross-sectional TEM imaging 
using a standard lift-out technique in a focused ion beam 
microscope [29]. Following a thinning step using 5 kV 
Ga+ ions, we performed a final cleaning and thinning 
using 900 eV Ar+ ions. 

An image of a two-angle Py film deposited on Gr/Cu, 
10 nm at 9 = 90° and 10 nm at 9 = 30°, is shown in 
Fig. 5. While the film is quite rough, it is continuous and, 
importantly, shows neither internal voids nor an apparent 
boundary between the two Py layers. This suggests that 


Aint (J/m^) 

Py/Cu, two-angle 

(6.6 ±0.6) X 10"^ 

Py/Gr/Cu, two-angle 

(5.6 ±0.6) X 10"^ 

Py/Cu, e = 30° 

(3.8 ±0.3) X 10“^ 

Py/Gr/Cu, 6 = 30° 

(3.5 ±0.3) X 10-* 


TABLE I. Anisotropy energies from FMR measurements. 


any holes in the first layer are filled in during the second 
deposition. The Py grains are fairly round and selected 
area diffraction patterns (not shown) indicate a modest 
(111) texture. 


IV. DISCUSSION 


Since our co-deposited Py/Cu and Py/Gr/Cu films 
have the same magnetization, similar perpendicular 
anisotropies, and nearly the same damping for films 
15 nm and thicker, the dramatic difference in damping 
for thinner films (see Fig. 1) must be attributed to the 
Gr interlayer. We do not expect other effects such as 
surface oxidation and surface roughness to be substan¬ 
tially different between the two types of samples. This 
leads to the surprising conclusion that a single layer of C 
atoms can block the transfer of spins between two met¬ 
als. Before offering an explanation for this effect, we first 
compare our results with other FMR measurements of 
Py films deposited on Gr. 

In [30] and [14], Gr grown by CVD on Cu foils was 
transferred to insulating or semiconducting substrates 
and Py was deposited by thermal evaporation to give 
Py/substrate and Py/Gr/substrate samples. Only a sin¬ 
gle Py thickness, 14 nm, was used for all samples in these 
two studies. In [30], for samples where the Py covered 
the entire Gr area, the Py damping was 1.9 times larger 
when Gr was present. This was interpreted as evidence 
for strong spin scattering in Gr itself, which would be con¬ 
trary to all expectations. However, Meff was also differ¬ 
ent for these samples (30% lower with Gr present) which 
suggests the difference in damping could be due to differ¬ 
ences in the Py itself rather than the presence of Gr. In 
[14], the same group made Py/Gr samples where the Gr 
extended « 10 pm past two edges of the Py film. Damp¬ 
ing for these sample was 1.1 times larger than for similar 
samples where the Py covered the entire Gr area. This 
small change in damping cannot be reliably attributed 
to spin pumping without further measurements to rule 
out changes in the Py. Unfortunately, no independent 
measurements of Ms or were reported in [30] or [14]. 
Further evidence of the need for caution comes from a 
study of 14 nm thick Co films deposited by evaporation 
onto Gr/Si02 and Si02 substrates [15]. Careful charac¬ 
terization by FMR, magnetometry, magnetic force mi¬ 
croscopy, and Kerr microscopy revealed such clear mag¬ 
netic differences between the samples that the authors, 
commendably, declined to draw any conclusions about 
spin pumping into Gr. 

In [13], Py and Pd were evaporated onto a Gr/Si02 
substrate and patterned into stripes separated by about 
1 pm of unpatterned Gr. When the Py was driven to 
process by an 9.62 GHz magnetic field, a voltage was de¬ 
tected in the Pd stripe that behaved as expected if a spin 
current was pumped from the Py into the Gr channel and 
then generated an inverse spin Hall voltage after diffusing 
to the Pd. FMR measurements at this same frequency 



5 


E 

c 


03 

(U 

03 


C 

(U 

CL 

■M 

c 

0 ) 

E 

o 


30 


20 


0 



0 10 20 30 40 50 

Nominal Py thickness c/q 


FIG. 3. Magnetic moment per unit area vs nominal Py thickness from VSM measurements. The gray shaded region indicates 
the initial 5 nm of Py deposited at glancing incidence, 6 = 90°. The rest of the Py was deposited a,t 9 — 30°. Horizontal error 
bars reflect an estimated 10% uncertainty in deposition rate, and the fits are weighted accordingly. 
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FIG. 4. Med vs. l/dpM with weighted linear fits for samples deposited at two angles and at a single angle {6 = 30°). 


showed the linewidth for Py/Gr/Si02 was 0.5 mT larger 
than for Py/Si02, for 25 nm thick Py. This difference 
was attributed to spin pumping from Py to Gr, but since 
AHq was not measured the difference could equally be 


due to differences in the Py itself. The fact that AHq for 
our 25 nm Py films differed by 0.9 mT, while a differed 
by only 1%, shows that FMR linewidth at a single fre¬ 
quency does not yield a reliable estimate for damping or 
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FIG. 5. TEM cross-section image of a Py/Gr/Cu sample with 
10 nm deposited a,t 0 — 90° and 10 nm deposited at 0 = 30°. 


spin pumping. 

The limitations of these previous experiments can be 
avoided by varying both FMR frequency and FM thick¬ 
ness over a wide range. The frequency range reveals what 
fraction of FMR linewidth is due to damping (a) as op¬ 
posed to inhomogeneity (AHq). The plot of a vs. (ipM 
allows a clear separation of the intrinsic and spin pump¬ 
ing contributions to a at each thickness. 



FIG. 6. Fit to damping vs. thickness for Py/Cu samples 
using Eq. 1. 

For our Py/Cu samples, the spin pumping expression 
in Eq. 1 can be fit to the data shown in Fig. 1 and re¬ 
produced in Fig. 6. For this ht we used the value of 
calculated from the data in Fig. 4. We also added an 
adjustable thickness representing a Py “dead layer” 
that does not contribute to spin pumping. Finally, we 
added an adjustable intrinsic damping oq to account for 
non-zero damping at large thickness. The fit then yields 
ao = 0.0074±0.0001,Geff = (8.6±1.6) x . ^-2^ 

and = 4. 7 ± 0.7 nm, and is shown by the dashed curve 
in Fig. 6. The uncertainty in is largely due to uncer¬ 
tainty in the thickness of the thinnest film. 


As described in Appendix A, we use a model for GeS 
that combines an interfacial spin mixing conductance G-i-j, 
and an external conductance Gext as resistors in series, 
i.e.. 


^ 

^eff ^ext 

The interfacial contribution is of order 10^® 17“^ • m~^ 
(see Appendix A), so Gefi in our case is dominated by 
Gext and insensitive to the value of G-f-j,. We can therefore 
neglect G^^ and, using the expression for Gext given in 
Appendix A and a bulk conductivity for Cu of cr = 5.8 x 
■ m“^, we infer a spin diffusion length in Cu of 
As = 310 ± 60 nm. This value is roughly consistent with 
the various room temperature values compiled in [31]. 

While we are unable to find prior measurements of spin 
pumping dead layers at Py/Cu interfaces, measurements 
on other FM/NM systems show much thinner dead lay¬ 
ers, ^ 1 nm [32]. This could be due to differences in 
the Py/Gr interface due to our unusual deposition con¬ 
ditions. Additionally, our 1-D spin pumping model does 
not account for lateral sample inhomogeneity, such as 
that from roughness observed in Fig. 5, but the effect of 
in-plane spin currents on our parameters is unknown. 

A likely explanation for the dramatic suppression of 
spin pumping by a single Gr layer is the electrical con¬ 
tact resistance between Gr and Gu. Measurements for 
several metals deposited on Gr gave a value of 1/Gc ~ 
10“® n • m^ for Cu and values of « 2 x 10“^^ fl ■ m^ for 
Ni, Ti, Pd, and Pt [16]. Adding this term to Eq. 4, we ex¬ 
pect Geff « Gc « 10® 11“^ ■ m“® for Py/Gr/Cu, which is 
nearly 5 orders of magnitude smaller than Ggff for Py/Cu. 
Thus the Gr/Cu interface presents a large barrier to spin 
flow, reducing the net flow of spin current away from the 
Py to the point where even our thinnest films show no 
enhanced damping. 

Even the lowest metal-Gr contact resistances are much 
larger than l/G-i-j., which raises the question of whether 
spin pumping into Gr is possible under any conditions. 
Although the values of and Ag inferred in [13] depend 
entirely on a questionable value for Aagp, the voltage de¬ 
tected at the Pd in those experiments is strong evidence 
that spins did flow from the Py to the Pd. Given this 
evidence, we believe it is likely that in our Py/Gr/Cu 
samples there is significant spin pumping into the Gr, 
but also a large spin accumulation because of the barrier 
at the Gr/Cu interface and therefore a large diffusion of 
spins back into the Py. A quantitative understanding of 
this situation, in which the same monolayer of C atoms 
participates in both the Py/Gr and the Gr/Cu interfaces, 
is a rich topic for further investigation. 


V. CONCLUSION 

In summary, we prepared Py/Cu and Py/Gr/Cu sam¬ 
ples with varying Py thickness dpM) using a two-angle 
sputtering method and carefully verifying that neither Gr 
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nor Py were degraded. We used FMR to measure damp¬ 
ing a for both series of samples. The Py/Cu samples 
showed an increase in a with decreasing dpM, as expected 
from spin pumping. Fitting the standard model of spin 
pumping showed GeS was dominated by Gext of the Cu 
rather than the interfacial conductance G-f-j,, which pre¬ 
cludes determining a value for the latter. The Py/Gr/Cu 
samples, despite having the same Py properties, showed 
no change in damping with thickness. This implies a 
strong suppression of spin flow away from the Py, which 
we attribute to the large contact resistance at the Gr/Cu 
interface. 

Our results are consistent with the evidence for spin 
pumping from Py into Gr reported in [13], where an in¬ 
verse spin Hall voltage was detected at a separate Pd elec¬ 
trode. However, we emphasize that quantitative analysis 
of spin pumping data requires independent characteriza¬ 
tion of magnetic properties and a consideration of the 
effect of all relevant conductances, not only G^j,. Our 
results are not consistent with the interpretation offered 
in [30] that Gr itself is a strong spin absorber. We ex¬ 
pect that magnetometry and/or a thickness series would 
show that differences in Py properties, rather than spin 
pumping, are responsible for the change in damping in 
this case. 


VI. METHODS 

A. Cu film deposition and Gr CVD 

Epitaxial, crystalline Gu(lll) thin films on sapphire 
were made as described in [33[. In brief, 500 nm of 
Gu was sputtered onto 50 mm wafers of q;-A 1203 ( 0001 ) 
held at 65 °C. During the annealing that preceded Gr 
growth (see below), secondary grain growth resulted in 
exclusively Cu(lll) grains > 2 cm across. These large 
grains prevent dewetting of the Gu film during Gr GVD 
at higher temperatures [33[. One GU/AI 2 O 3 wafer was 
removed after annealing only and used for Py deposition 
directly on Gu, while a second wafer continued through 
the Gr growth step. 

Gr GVD was performed in a hot-wall quartz tube fur¬ 
nace with a diameter of 76 mm. After evacuating the 
tube with a dry pump to a pressure of 0.67Pa (5 mTorr), 
the rest of the process was performed at a total pressure 
of SOkPa (600 Torr). The wafer was heated in a flow 
of 7 seem H 2 and 2100 seem Ar (H 2 partial pressure of 
0.3kPa (2 Torr)), initially ramping to 920 °C in about 
15 min and then moving to 1000 °G at 2°C/min to al¬ 
low complete Gu grain growth. The temperature was 
then increased to 1060°G and the gas flows were changed 
to 58 seem H 2 and 3600 seem Ar (H 2 partial pressure 
of 1.3kPa (9.5 Torr)). After 10 min under these condi¬ 
tions, Gr growth was initiated by adding 18 seem of 0.2% 
GH 4 in Ar (CH 4 partial pressure of 0.8 Pa (6 mTorr)). 
Growth conditions were maintained for 3 h before rapid 
cooling by opening the furnace lid. The flow of GH 4 was 


stopped below 400° G and the flow of H 2 was stopped 
below 100 °C. 

The resulting Gr hlms were characterized by optical 
microscopy to check for completeness of growth and Ra¬ 
man spectroscopy to confirm their quality. For shorter 
growth times, for which Gr coverage was incomplete, 
we found this growth recipe consistently gave compact 
hexagonal Gr growth domains with nucleation sites sep¬ 
arated by « 100 pm. 

Both wafers were cut into 6 mm x 8 mm chips using 
a diamond scribe. 


B. Py sputter deposition 

For each thickness, Py was deposited simultaneously 
onto bare Gu and Gr/Cu samples. To remove Cu oxide 
and adventitious contamination, samples were soaked in 
glacial acetic acid for 1 min and rinsed in deionized wa¬ 
ter just prior to being placed in the load lock of a de¬ 
position chamber with a base pressure w 4 x 10“® Pa 
(3 X 10“®Torr). The samples were heated to « 200 °C 
to remove residual water and then cooled in vacuum to 
« 50 °C before deposition. Depositions at both 6 = 90° 
and 6 = 30° were performed with the source running 
at 200 W and w 430 V, in 1.1 mTorr of Ar gas, and 
with a target-to-sample distance of 23 cm. The deposi¬ 
tion rates, measured using atomic force microscopy for 
« 10 nm thick hlms, were 0.03 nm/s for 6 = 90° and 
0.09nm/s for 6 = 30°. 

C. Magnetometry 

The VSM measurements were performed at room tem¬ 
perature with the sample mounted in a plastic straw, us¬ 
ing a frequency of 40 Hz, and an oscillation amplitude of 
2 mm. The magnetic held, applied in the plane of the 
sample, was swept over a range of hqH = ±4 T to allow 
a linear ht to the diamagnetic background, due primarily 
to the AI 2 O 3 substrate. The moment of the Py layer was 
taken as the H = 0 intercept value of a linear ht from 
1 T to 4 T. 


D. Raman spectroscopy 

We performed Raman spectroscopy in a home-built 
system using a 532 nm diode laser with a power of 
« 3.5 mW and a spot size < 10/rm. Peak htting was 
done after subtracting a smooth ht to the surrounding 
background from Cu huorescence. 
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APPENDIX A: MODEL FOR SPIN PUMPING IN 
MULTI-LAYER STRUCTURES 

Following [32], we model the spin transport of a FM 
layer in contact with one or more NM layers as an electri¬ 
cal circuit of resistors in series. Fig. 7 shows the simplest 
case of one interfacial resistance, l/G-i-j,, and one exter¬ 
nal resistance, 1/Gext, The effective conductance for spin 
current flowing from the FM to the spin sink is simply 

1 _ I 1 

Since these conductances describe spin transport perpen¬ 
dicular to the FM/NM interface, they are defined per 
unit area of the interface and have units of -m”^. 
In the full theory of spin pumping [1], G-i-j, is a com¬ 
plex quantity whose real part affects damping and whose 
imaginary part affects the FM precession frequency. For 
conventional FM metals such as Ni, Fe, and Co, the ef¬ 
fect on frequency is negligible [34] and G-j-j, is commonly 
used to refer only to the real part. 

1/G„ 

FM Spin 

si nk 

FIG. 7. Series resistor model for spin transport between a 
FM and a NM spin sink. 

The interfacial conductance G^j,, often called the “spin¬ 
mixing” conductance, determines how readily spins can 
move across the FM/NM interface. This process is gov¬ 
erned by the exchange interaction, which operates over 
a length scale of a few nanometers. For metallic in¬ 
terfaces, the theoretical upper bound on G^j. [7] is the 
Sharvin conductance [35] (for a single spin direction and 


per unit area), corresponding to perfect transmission for 
every conductance channel at the interface, 


Gs 


h Jtt ’ 


where fcp is the Fermi wavevector. Taking Gu as a typ¬ 
ical metal, with = 1.36 x 10^°m“^, we find Gs « 
0.5 X 10^® ■ m“^. Both experimental and theoretical 

results for Py/Cu interfaces (and several other FM/NM 
combinations) are consistent with this limit [5-9], with 
G-f-j, in the range of (0.5 ± 0.1) x 10^® • m“^. 

The external conductance Gext determines how readily 
spins flow away from the interface once they have entered 
the NM material. For a single NM layer of thickness (Inm, 
bulk electrical conductivity a, and spin diffusion length 
As, Gext is given in [32] as 

1 O' (Inm 

LText — ^— t&llil —; . 

CA. r) \ \ 

Z Ag Ag 


For cInm ^ 2As, which is often the case, the external 
spin conductance is simply determined by the ratio of 
conductivity and spin diffusion length, Gext ~ o’/ 2 As. 

In the limit of Gext ^ when spins diffuse rapidly 
away from the interface, Geg G-i-j, and a measurement 
of Aofsp can be used to directly infer a value of G^^ from 
Eq. 1. This limit applies to NM materials with large a 
and/or small As, such as Pt and Pd [32]. However, other 
materials such as Ta and Gu [32] are in the opposite limit, 
where Geg ~ Gext and Aasp is not sensitive to G^^. 


APPENDIX B: FULL FMR RESULTS 

Table II shows all parameters obtained by fitting the 
FMR data for each sample. These results are plotted in 
Figs. 8 , 9, and 10. Values of g come from fits to Equation 
2 , where 7 = 


0.030 

0.025 

t 0.020 

O 

^ 0.015 

O 

^ 0.010 

0.005 

0.000 


FIG. 8 . Plot of AHo vs Py thickness. 
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do (nm) 

goMeff (T) 

9 

fioAHo (mT) 

a 

Cu-tGr 

5-t3 

0.6879 

2.089 

31.6 

0.0073+0.0001 


5+3 

0.7918 

2.098 

18.2 

0.0077+0.0002 


5+10 

0.9002 

2.108 

11.9 

0.0072+0.0005 


5+20 

0.9509 

2.097 

1.7 

0.0084+0.0007 


5+40 

0.9772 

2.102 

3.2 

0.0089+0.0005 

Gu 

5+3 

0.7298 

2.098 

13.4 

0.0138+0.0003 


5+5 

0.7654 

2.099 

4.9 

0.0089+0.0001 


5+10 

0.8807 

2.102 

0.9 

0.0078+0.0001 


5+20 

0.9381 

2.099 

0.8 

0.0075+0.0002 


5+40 

0.9725 

2.104 

1.1 

0.0076+0.0001 


TABLE 11. Parameters from FMR fits of Ji'res(/) and AH{f). 
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FIG. 9. Plot of Meff vs Py thickness. 
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APPENDIX C: GRAPHENE DAMAGE 
ESTIMATES 

Expected Damage Rate from Ar Neutrals 

The simulation results in [20] can be applied to our 
sputtering conditions to estimate the flux of Ar neutrals 
arriving at the substrate with sufficient energy to damage 
the Gr. In particular, the results in Fig. 7 of [20] for 400 
eV Ar+ ions hitting a Cu target are a good estimate for 
our case of 430 eV Ar+ hitting a Ni 8 oFe 2 o target. The 
probability of an Ar neutral reflected normal to the target 
and with sufficient energy to damage Gr is approximately 
2 X 10“"^ per unit solid angle for each Ar+ ion that hits 
the target. Our sputtering source operates at a current of 
0.47 A, corresponding to 2.9 x 10^® Ar+ ions per second. 
Our 6 mm x 8 mm substrate, oriented at 0 = 30° and 
at a distance of 230 mm from the target, subtends a 
solid angle of 7.9 x 10“"'^ sr. The result in Fig. 21 of [20] 
shows that there is negligible scattering of Ar neutrals 
for our pressure-distance product of 35 Pa • mm. This 
yields a flux of damaging Ar neutrals at our substrate 
of about 0.009 nm“^ • s“^, remarkably close to the value 
of 0.005 nm“^ • s“^ inferred from the Raman spectrum of 
the 2 s exposure described in the main text. 


Raman Measurement of Gr Damage vs. Angle 

In order to evaluate Gr damage as a function of sput¬ 
tering angle, we deposited « 10 nm of Py at several angles 
between 9 — 30° and 9 = 90° and made Raman measure¬ 
ments. The results for Id/Iq and I^d/Iq are shown in 
Fig. 11. 


FIG. 10. Plot of g vs Py thickness. 
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Py deposition angie 0 



Py deposition angie 0 

FIG. 11. Raman peak ratios Iu/1g and ho/Ia vs Py angle 
of incidence 9. 
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